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SUMMAHT 



to analysis of a solid propellant ramjet which carries the 
propellant in the combustion chamber was made, and a procedure 
outlined for calculating the required burning rate of such a ramjet 
which develops constant thrust. The two factors which influenced 
this development were: 

1) The solid propellant ramjet's combustion chamber varies 
in si se with tine due to burning of the fuel. 

2) The 3olid propellant ramjet develops constant thrust. 

These factors wore coordinated with the analysis of the internal 
flow system of the ramjet, and an oppression for the required 
burning rate was derived. 
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«= area, square feet 

t= velocity of sound, feet per second 

« drag coefficient of combustion chamber 

» drag coefficient for incompressible flow 

*= 3 kin friction coefficient 

= specific heat at constant pressure 

(for air 0.24 lVHj/lb/°F) 

«= throat diameter, inches 

= internal diameter between individual grains, or charges, 
making up combustion surface. Inches 

*= drag force due to friction in duct 

= incompressible drag force due to friction in duct 

= system parameter, S * F/»p 

*3 fuel - air ratio 

= thrust, pounds 

*= acceleration of gravity (32.3 ft/sec ) 

» heating value of fuel, BBJ per pound 
** h9at added, 33SJ per pound 

= specific impulse, seconds 

jb length of combustion chamber, inches 

» mass flow rate of air, pounds per second 

=: fuel mss flow rate, pounds per second 

* naeh number (v/a) 



number of Individual grains making up combustion 
surface 

static pressure, pounds par square inch absolute 

dynamic pressure, pounds per square inch absolute 

gas constant, feet per °F (for air 53.3) 

specific fbel consumption, weight of fuel consumed 
per hour per pound of thrust 

Idoal specific fuel consumption, weight of fuel 
consumed per hour per pound of thrust 

time of burning, seconds 

temperature, °F absolute 

total temperature, °? absolute 

velocity within duct, feet per second 

weight of fuel consumed, pounds 

ratio of specific heats (for air 1.40) 

density, pounds per cubic foot 

stations in internal flow system shown in Figure 1 
fuel 



throat 



i. nwnonocTios 



T!i 3 U 39 of solid fuel in a ramjet introduces a new problem. 

As the fuel bums, the effective combustion chamber area increases. 
This would normally result in a change in thrust; however, if the 
condition of constant thrust is deposed, thsn this change in 
combustion chamber area would vary the amount of heat required. 

The change in heat required nay be expressed in terms of the burning 
rate required for a solid propellent ramjet developing constant 
thrust. 

References 1 and 2, which are representative of the numerous 
analyses in this field, were concerned with the internal flow through 
a ramjet having a constant combustion chamber area, and no restriction 
on the thrust. These analyses were developed for subsonic conditions, 
and a United heat addition. 

In the present investigation, a theoretical analysis was made 
for supersonic conditions of a solid fuel ramjet, which carries the 
propellant in the combustion chamber, developing constant thrust. 

Thi3 latter requirement fixes the exit conditions and permits 
calculation of the internal flow from the rear of the ramjet to the 
combustion chamber inlet. In this manner, the pressure loss 
across the combustion chamber and nozzle was computed. The pressure 
rise across the diffuser may be determined from the hnown inlet 
conditions. 3y equating the pressure rise to the pressure loss, 
unique determination of the heat required for constant thrust can be 
made. The amount of heat developed is a function of the fuel and 



its distribution, and from thi3, tho required burning rato nay 
bo de to mined. 

I'.-csults of this -.U'iysis '*'(? pros-nted in Trvt II of the 
invent! j tion in the form of v ■orbing equations, some of which have 
boo:: solved ?nd presented on ch- rts. Theoe equations and charts 
have boon consolidated, and the required bumirr; rate or 3 i>ocific 
fuel consumption is presented in Figure 10. Part III illustrates 
tho method of solving a numerical problem. 
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Figure 1 



Schematic of a solid propellant ramjet which 
carries its fhel in the combustion chamber. 



II. ANALYSIS 07 IH22SML HOW SYSTUiS 



51o u characteristics of the internal flow system cast bo 
calculated go that the state of air entering the combustion chamber, 
the flow change across the combustion chanher, the net drag of the 
internal flow, and the required geometry of the ramjet can he 
determined. In determining the internal flow fron tho roar end of 
the ramjet, the follord-ng relationships mst ho established; 

(A) nozzle throat condition to ambient condition; (B) flow changes 
across the nozzle; and (c) flow changes across the combustion 
chamber. Those expressions nust then be related to (D) flow changos 
across the diffuser. 

One dimensional flow, unlfom velocity distribution across the 
duct, constant gas properties, and ideal expansion were assumed. 

Initial quantities assumed to bo known are as follows; 

(1) The desired thrust at designed altitude. 

(3) The heating value and density of the fuel. 

(3) Pesigned Mach number (supersonic) and atmospheric 

conditions p 0 , f> 0 , T 0 . 

(4) Skin friction coefficient. 

(5) Area of the throat of the nozzle. 

(6) Ratio of length of chamber to throat diameter. 

(7) The number of grains (N ) of fUel. 



A conplote investigatio \ sTiOuld include consideration of: 
l) lsentropic diffusion, and 2) nomal shock wave; and whether 
the nozzle Is sonic or not at the throat. This analysis outlines 
in detail the case of isentropic diffusion with a sonic nozzle. 

A. Slow Relation Sotwosn Hosslo 3Mt rad Aablont Conditions 

.to e;caot relation between no sale exit condition and that of 
the nnbient surroundings nay bo obtained fren the expression for 
thrust. The thrust of a ran jet io givon by the increase in 
nonantun of tho internal flow, or 

F = is ( v 4 - v Q ). 

The mss rate of flow at the nozzle throat is 

a “ p\ v t H** ft a t A t since M*«l 
The thrust nay then be expressed as 



S ' = ft a t A t ( V 4- V 0 J 







.All of tho factors on tho loft hand 3ide of the equation aro 
Inoun initially, therefore a new pareneter, E , is introduced 
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Then by tho perfect gas law 
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Using the isentropic relation 
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Solving for u } 
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( 1 ) 



(1A) 



is not a known parsnoter, hut it can he related to the tomperaturo 
rise. The heat liberated is equal to the total tenrperature rise, or 



( 2 ) 
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Tho conprossion from freo strean to station 2 at the conbustion 
chamber inlot i3 an isentropic one, and the state 2 is assunod to 
be a stagnation state, where 



Similarly the expansion from state 3 to state 4 is an isentropic one 
and 

lit - 



Equation (2) becones 

J±c _ j 
7l t ~ Cp lot: 



(3) 



Using the isentropic relation 



/4<- - T 4 C I + 3. ) 



equation (3) becones 
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(4) 



Squating (1A) and (4), and solving for 
following expression. 
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, results in the 
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This equation con 
v o solved In tho form 
of a graph and would bo 
oinil,. .r to that sho-.jn 
In rigor 3 2. 




The Ideal valuo of heat requlrod ( a ramjet with no 

losses can ho obtained by entei'ing Figure 2 with H 4 « H 0 » Aa losses 

develop, the exit mach rronber % decreases, and a larger value of 
A H 

— =— is required to maintain constant thrust. 

Cp i.r 



B. Flow Changes Across Tho ITosslo 

In development of the analysis from the rear of the ramjet, tho 
flow across the nozzle must be expressed in terns of the combustion 
chamber outlet. The pressure ratio across tho nozzle is 



Jr 




or, in terns of mach numbers, 

p 3 __ < r 

lO V I +- Ml I 



hence 

Ti I -t- 

T* I + ^ ^ 



(?) 



Ifcr the continuity equation 
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or 
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Mao, 




Solving for t^. , 
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Substitution of this relation into equation (6) yields 
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(9) 



and Hg any he obtained hy equating (7) 
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A relation between 
and (8) 



(rT J>c. 



I + 



( 10 ) 



In o-ia’ticn (9), the torn frT I><- is a known quantity and 

haa physic?! significance, whereas Mg is an unknown value and is 

not generally used as a parameter. Shis analysis will thorofore 

d 

so developed in toms of ^ ; however, because of oase of 

conversion, as seen in equation (10), solutions to all equations 

wore derived in toms of Mg. The final result is expressed in the 

d. _ 

known paranotor, . 

Equation (9) nay be written as 



2s mow Changes Across The Combustion Ghgnber 

The flow through the conbustion chamber is considered nest. The 
c)\annel or duct has a constant cross section in which friction losses 
occur. As the fuel bums uniformly, the cross section of the duct 
increases, and a change in heat added is required to maintain constant 
thrust. In keeping with the procedure outlined, all relations will 
be stated in terms of conbustion chariber outlet conditions. 

The pressure drop between stations 2 and 3 (at conbustion 
chamber inlet and outlet) is obtained by applying the nonentoE relation 



x 




( 11 ) 



befcwoon these stations: 



where % is the drag forco duo to friction within tho combustion 
chamber, end Aq = Ag. 

From equation (12) it is evident that the pressure drop la 
composed of a component associated with tho frictional drag force, 
and a corrponent duo to the momentum increase in the duct. is 
pointed out in Reference (l), the former component is frequently 
neglected hut becomes very important when the velocity in the 
duct is large. Tho drag term is presented in the fom of a drag 
coefficient, Cj>f , where 

d>f 

Cj> f = (13) 



If the incompressible flow is taken as a basis for the unheated low 
speed condition, the drag coefficient becomes 



Cr>u - 



(- 



3>* 

<3*. Av 



) 



\ 



(14) 



end is determined by the geometry of tho ramjet sad the Reynolds number. 
There is nn increase in the dynamic pressure when heat is applied , 
resulting in en Increase in the incompressible drag coefficient. 

If the dynamic pressure varies linearly with the length of the 
combustion chamber, an average value may be assumed, and equation (14) 
becomes 



If the incompressible drag is now expressed in terms of the skin 
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frlction coefficient (a known value), C Decodes 



or 



hence 



C 



J>fA 



, r tv/* 

Cf T[ J — J Lrvrr £>c 

J_ /V TT £>c 



4. 



C Df a - 4 - j> 0 W 



/* / f (iVi t 7 /i 

2> fi - 4 it c f r|_ — 2 —]^ 



(15) 



Since the heating effect has been accounted, for in the basic drag 
relation, the two drag equations (13) and (15) nay be equated 

Cm Cl f% /), =fjc f (<£ + |)4, 



or 



Cm = ( l+ %) 



(16) 



With this definition of the drag coefficient, equation (12) beconea 



p - + p : vC - p 3 + ^ v , 1 + ^ X c f ^ (<-*- /fy ) 



or 



- I + 



X c f ('-"%) + 0-/y (V )] 
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or, in terns of tho dienoter ratio, 3 inco D c varies, 




* - 6 
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( 17 ) 



In order to evaluate the pressure ratio, tho density ratio 
must first he obtained. This ratio nay he determined by simultaneous 
solution of the energy equation, 

~t + Cp lx t 4H - ^- + C P J 3 



the continuity equation, 

p 1 **. ^ C* V 2 

the perfect gas relation, 

p = £ R.T 

and the momentum relation, equation (12). The final result is 
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Solution of equation (17) provides 
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for, 
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f 0 *- _ _ Px. 

Po " Po p, (19) 



Thus the pressure drop across the combustion chamber and nozzle has 
been calculated, and nu3t equal the pressure ri3e across the 
diffuser. 

Solution of equation (17) also permit s determination of f 

for by tho continuity equation, 

(>x ax M*. = a 3 Ms 



or 



Also, 



&)* = %) 
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Combining the last two equations, and simplifying, 
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( 20 ) 



The pressure rise across the diffusor as computed from the rear 
of the ramjet was plotted vs the combustion chamber inlet nach 
number as calculated in equation (20), (see Figure 3). 
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A plot of combustion chamber inlet naeh number and heat 
required for constant thrust as developed in equations (18) 

and (30) is shown in Figure 4. 




Figure 4 
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2z Flow Changes Across Tho Diffuser 

In an analysis of a ramjet haring a designed suporsonic mch 
number, there are two types of diffusion to be considered - namely, 
i3entropic and normal shock wave diffusion. In this investigation, 
diffusion was considered to be isentropic, and the diffuser 100$ 
efficient. Tho pressure ri3o across the diffuser is 




( 21 ) 



how, the pressure rise across the diffuser as computed from the 
front end of the ramjet was plotted in Figure 5. 

N*. vs. A- for constant E, Mq, 



r\ 

Po 

Figure 5 

p 

Unique determination of W*. and for constant 

thrust was obtained by superimposing Figure 5 on Figure 3 and 
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noting the points of intersection nn shown in Figure 

vs - for constant E, 1^ 
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de t ermination of Iteming Rate 

The condition of constant thrust has been satisfied by the 
simultaneous solution of equations (19), (20), and (21). With the 
valuo of Mi obtained at the intersection of the two curves as 
a parameter, the value of heat required ) was determined 



from Figure 4 for a given value of 



After obtaining values 



^ H d 

of for various fn -be. ♦ a chart was nest constructed of 

cL 



(w be 




|7TL 



Figure 7 
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7- is chart was used later to obtain the heat required for an 
elemental volume. In considering m elemental volume, the 
weight of fuel consumed and the tine of consumption must bo 
calculated. 

To accomplish this, the mass rate of flew is written as 

/V-iT^l fl a*. 

m = (l a, Mx/U = fi a. A?*. ~ * Trsr 



where 



Then, 



^3. 

^ a « 









H I 




(22) 



The fuel-air ratio say be determined from 
= K 

P = m 

The H was obtained from Figure 6, and h is a given value, 
therefore 

m P = m f 

The final relation between the weight of fuel burned, the tine 
of burning, and the changes in diameter may be obtained by 
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concidcrlng a small elemental volume of fuel, The volume ( V ) 
■at tine "t is 

V, = L Nv i><_ ( 

where A D c is the change in fuel diameter. At tine t + it , 



V, 



i> t Ct + 4t)= .be (.t) t Abe. - 2>c (*) + L/VTT hc(0 



(24) 
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Or, in toms of ~~j— , 
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Hence, solving for V, f 

r 1 ^ i 1 

V, - LNTTJ>t(t)^ $ 6c K) - )J 



(t >4t) 



V, = U L At V 



Ac ( *)J 
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(25) 



The weight of the fuel consumed, w p , nay be obtained fron 



Wp — V, 
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vho tiio of burning of this elcnentnl volune Is 

- m F 

A plot of the weight of fuel consumed V3. disaster ratio and 
vs. tine of burning was nade in Figures 8 and 9. 





Figure 9 



So show the curve in detail, the specific fuel consuaption 
ninu3 the ideal specific fuel consunption was plotted. Bio specific 
fuel consumption is 



3feoo f 3<*oo v»\p 

S = r = f 



(26) 



For the ideal case in which there are no losses, 



v-p 0 Ar 



TT + t 
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I i- 
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(27) 



Tho v ine of CpToc 5‘ obtaijwd, from Figure 2 entering with 
M+ ** A'Iq ?.a a parrueter. Then 

3 boo a H m g. 

S 0 = K F 



where m is as expressed in equation (27). 

A plot was then constructed of ( S - So) vs. tine of burning 
(see Figure 10). 
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Figure 10 
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III . A . U.l 1 'JO. 0. i nOAL J ; 

1:1 a numerical problem to determine the burning rate of a 
solid fuol ran jot (which carries its fuol in tho combustion 
chonbor) developing constant thrust, tho following data. is given: 



n o = 


2 


ll = 


14,230.65 eru//fe. 


H “ 25 


3 2= 


2 


= 


140.5 1B /ft 3 


0 f = 0.003 


If := 


1.4 


At « 


1 ft 2 


¥ - 60 


P 0 = 


4.362 psia. 


y 

o 


411 .7° ft 


F = 1758 lb. 



Squat ion (5) has been solved for those values, and is shorn in 

Figure 11. From this chart, the ideal value of heat required for a 

ran jet with no losses is 1.019. As the losses develop, the exit 

A H 

nach number is decreased and a larger value of 7 ~=" is 
required. 

d- 

Hie relation between fc and as stressed in 

equation (10 ) i9 shown in Figure 12. 

Solutions of equations (19) and (20) for various values of 

cL 

are shown in Figures 13A, 13B, and 130. It is to be recalled that 
those solutions were calculated from the rear of the ramjet. 
Equation (2l), calculated fron tho front end, i3 plotted on these 
sane figures, and the intersection of these two curves uniquely 
determine — and Mi- for a constant thrust, solid 

Po 






propellant ramjet 
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The heat required to maintain this constant thrust was 
plotted vs. tho combustion chamber inlot much number, and is 
shown in Figures 14A., 14JJ, and 140. Using the uniquely determined 
valtte of Mi. as tho entering parameter, the heat required was 
obtained as a function of the diameter ratio ) , and is 

shown in Figure 15. 

Using the process of an elemental volume as described in 
Section T. t and solving equation (25), a relation between the wei^t 
of the fuel consumed and the diameter ratio was determined. This 
relatlon^is shown in Figure 16. She weight of fuel consumed was 
plotted vs. the burning tine in Figure 17. The slope of this curve 
is the required burning rate. Unfortunately, for the values chosen, 
the burning rate i3 relatively constant. Actually, the required 
burning rate should decrease with an increase in burning tine. So 
show tho relation between tho weight of fuel consumed and the 
burning time in greater detail, tho specific fuel consumption 
was calculated from equation (26) and the value of (S - S Q ) was 
plotted vs. tine of burning and is shown in Figure 18. This curve 
shows a change of approximately fifty percent in (S - Sq) over tire 
range of burning tine. 
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